Towards Phonon-Like Excitations of Instanton Liquid 
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The model of (anti)-instanton liquid correctly seizes many nonperturbative phenomena and 
important vacuum features such as chiral symmetry breaking, the presence of gluon condensate 
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Abstract 

The phonon-like excitations of (anti)-instanton {II) liquid due to adiabatic variations 
■ of vacuum wave functions are studied in this paper. The kinetic energy term is found and 

Q\ ! the proper effective Lagrangian for such excitations is evaluated. The properties of their 

spectrum, while corresponding masses are defined by Aqqd, are investigated. 

£ : 

S: 
<n : 

(N . 

> '. 

and topological susceptibility It is usually supposed that the corresponding functional 

integral in this approach is saturated by quasi-classical configurations close to the exact solu- 
tions of the Yang- Mills equations (the Euclidean solutions called the (anti)-instantons) and the 
wave function of vacuum, being homogeneous in metric space, is properly reproduced by aver- 
aging over their collective coordinates. In the // liquid approach one takes the superposition 
ansatz of the pseudo-particle (PP) fields as one of the simplest relevant approximations to the 
Qh| 'genuine' vacuum configuration 
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Here A fl (x; , ~fi) denotes the field of a singled (anti)-instanton in singular gauge with 4N C (for 
the SU(N C ) group) coordinates 7 = (p, z, Q), of size p with the coordinate of its centre z, Q 
H ! as its colour orientation and 

Al(x- 1 ) = -W b f ]b »v y - 2 -^— 2 , y = x-z, (2) 

where 7] is the 't Hooft symbol 0. For anti- instanton fj — > 77 (making the choice of the singular 
gauge allows us to sum up the solution preserving the asymptotic behaviour). For simplicity 
we shall not introduce different symbols for instanton and anti-instanton, and then in the 
superposition of Eq.(|l]) N implies the PP total number in the 4- volume V system with the 
density n = N/V. The action of the instanton liquid model is introduced by the following 
functional 

(S) = [ d 4 x [ dp n(p)s(p) . (3) 



The integration should be performed over the system volume along with averaging the one 
instanton action s(p) weighted by instanton size distribution function n(p). Then an action per 
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one instanton is given by well-known expression 

sx(p) = M + 5 ln(Ap) - ln£ 2JV ° + / d Pl n{p x )p\ , (4) 

with the Gell-Mann-Low beta function ft(p) = — In Cn c — 61n(Ap) , A = A-^=r = 0.92Apy. , and 
constant Cn c depending on the regularization scheme, here Cn c ~ ^In^-i^In 8 ^)} ' ^ = T^ c ' 
and the parameters ft = ft{p) and ft = ft + In Cn c are the ft function values at the fixed quantity 
of p (average instanton size). 

Some terms (inperfect p dependence) of Eq. (|J) could be obtained in the classical Yang-Mills 
theory with one-loop (quantum) corrections taken into account and resulting in a modification 
of coupling constant g at the distinct scales. Indeed, the first term is the one instanton ac- 
tion 8ir 2 /g 2 with the p-dependence of g corrected. The last term of Eq.(||) describes the pair 
interaction of PP's in the instanton ensemble with the constant £ characterizing, in a sense, 
the intensity of interaction £ 2 = ^ Jt'Li ^ 2 - The smallness of characteristic instanton liquid 
parameter rap 4 ('packing fraction') allows us to drop out the p dependence of the /3-function. 
The logarithmic terms correspond to the pre-exponential factor contribution to the functional 
integral and are of a genuine quantum nature. In the second term the p 5 factor makes the 
integration elements dp and d^z dimensionless. Finally, the third term is a square root of the 
one instanton action ft raised to the power 4iV c . The latter is just the zero mode number of the 
one instanton solution and the corresponding p dependence may be again omitted because of 
small logarifmic contribution. 

Taking the exponential form for the distribution function over the instanton action n(p) ~ 
e -si(p) ; we obtain directly from Eq.(|3|) the self-consistent description of equilibrium state of 
instanton liquid with the well-known ground state 

u(p) = p-t^e-M-"*/? , v = 1(6-4) = , n = J dp n(p) , J = J dp p 2 n{p)/n. 

The distribution u(p) has obvious physical meaning, namely, the quantity d 4 x dp p(p) is 
proportional to the probability to find an instanton of size p at some point of a volume element 
d A x. At small p the behaviour of distribution function is stipulated by the quantum-mechanical 
uncertainty principle preventing a solution being compressed at a point (radiative correction). 
At large p the constraint comes from the repulsive interaction between the PP's which is 
amplified with (anti)-instanton size increasing. 

Deriving Eq.(|3|) we should average over the instanton positions in a metric space. It is clear 
that the characteristic size, which has to be taken into account, should be larger enough than 
the mean instanton size p. But it should not be too large because the far ranged fragments 
of instanton liquid are not causally dependent. The vacuum wave function is expected to 
be homogeneous on this scale L > R (R is an average distance between the PP's). Let us 
remind that each PP contributes to the functional integral with the weight factor proportional 
to ~ 1/V, V = L A . The characteristic configuration saturating the functional integral is taken 
as the superposition ([[]) with N pseudoparticles in the volume V. If one supposes that the PP 

4 This argument corresponds to the maximum principle of Q]. Approaching the functional (|[) as a local 
form (S) = J dp si(p)n(p)/n where si(p) = f3(p) + 51n(Ap) - ln/3 2Arc + /3^ 2 p 2 np2, using n(p) = Ce~ s ^ with 
the constant C as a distribution function (actually it makes the problem self-consistent because an equillibrium 
distribution function should be dependent on an action only) and taking the variation of (S) — (Si) — J dp {s(p) — 
s i(p)} e ^ p V n ov er s(p) one may come to the result s(p) = Si(p) + const keeping into the mind an arbitrary 
normalization. 
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number in an ensemble is still appropriate to consider them separately, then denoting AN(pi) 
as the PP number of size p G (pi, pi + Ap), K as the number of partitions within the interval 
(pi, pf), Eq.(|l]) may be rewritten in the following form 

K AN(pi) 

A ^ X ) = Y. £ An(x;i,Jj) , (5) 
i=i j=i 

where A^x; i, jj) is the (anti)-instanton solution with the calibrated size and 7 = (z, Q) stands 
for the coordinate of its centre and colour orientation. By definition Y^f=i AA^(pj) = N. Further, 
introducing the distribution function n{p) = y, and normalizing it as Y^f=\ n {Pi) Ap V = 

N (in the continual limit Ap — > it is valid V J dp n(p) = N) one can calculate the classical 
action S c — | / d 4 x G 2 of this configuration averaging over the instanton positions in the 
metric and colour spaces. As a result (with the superposition ansatz Eq. ([[])) one instanton 
actions and the PP pair interactions only contribute to the average system action 



(Sc) = fi / S <=J d * z J d P <P) {^f + / d Pl n(Pi)pl] , (6) 



where dfl is a measure in the colour space with the unit normalization. As above mentioned 
(S c ) including one loop corrections will then contribute to the functional integral. 

It is easy to understand that Eq.@ describes properly even non-equilibrium states of the 
instanton liquid when the distribution function n(p) does not coincide with the ground state 
one p(p). Moreover, it allows us to generalize Eq.@ for the non-homogeneous liquid, when the 
size of the non-homogeneity obeys the obvious constraint A > L > p. 

In what follows, we study the excitations of II liquid induced by adiabatic dilatational 
deformations of the instanton solutions. Then, as the configurations saturating the functional 
integral we consider not the instanton solution itself but the quasizero modes which are para- 
metrically very close in the functional space (a direction does exist where the action varies 
slowly) to the zero modes. The guiding idea of selecting a deformation originates from trans- 
parent observation. The deformations measured in units of the action (here g, p are the 
generalized coordinate and momentum) have a physical meaning only. However, the instantons 
are characterized by 'static' coordinates 7 and, therefore, need to appoint the conjugated mo- 
menta. It looks quite natural for the variable, for example, p to introduce those as p = dp/dx^. 

Let us calculate first of all the corrections for the one-instanton action. Dealing with super- 
position ansatz Eq.(|l]) again, one should include additional contribution to the chromoelectric 
field 

G% = + 9°"nv ■ ( 7 ) 

with the first term of strength tensor (s.t.) corresponding to the contribution generated by the 
instanton profile 



~ 9 {y 2 + P 2 ) 2 I2 M p y 2 r]ailp y 2 ) ' 

and in adiabatic approximation the corrections have the form 

d 4 y p 

9a ~ — ^ P = ~Vaiv , 9 ^ttt p, g^j = , <7 i4 = — <7 4i , — 1, 2, 3 
dp g [y I + p 2 Y 
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Here we have justifiably ignored the terms of order O(p), 0(p 2 ). The adiabatic constraint 
9% ^ means that the variation of instanton size is much smaller than the characteristic 
transformation scale of the PP field, p <C 0(1). Then calculating the corrections for the action, 
it is reasonable to take out p beyond the integral and the one instanton contribution to the 
action turns out to be 

1 r 8tt 2 



I J d*xG' 2 u ^ — + C p + ^- p\ (8) 

where p should be taken as the mean rate of slow solution deformation at a characteristic 
instanton lifetime ~ p. For simplicity, one may take it in the centre of the instanton p(z). The 
constant C = (because the first variation of the action SS/ 5 A for the solution itself equals to 
zero). For the 'kinematical' K-term we have 

12tt 2 

« a . t . = —5- ■ 9 

g 2 

The overt p dependence of k is lacking because of the scale invariance. It arises with the 
renormalization of the coupling constant (in a regular gauge the result is the same). Being 
within the ansatz (H) we have considered only the corrections induced by the variation of 
strength tensors, but not those resulting from a possible variation of fields @. Bearing in mind 
the form of potentials in regular (r.g.) and singular (s.g.) gauges 

A« = ^ Vatlu d u My 2 + p 2 ) (r.g.); A a ^ d u \n (l + Q (s.g.), (10) 

we find the adiabatic corrections A^ = A^ + follows: 

al = -Va^ 2 P . o P ( r -S-) • (11) 
" g y 2 + p 2 

The substitution r] — ► —fj brings about the transition from regular gauge to a singular one. 
Using the admixture gr* of chromoelectric and chromomagnetic fields generated by the to 
saturate the functional integral as in Eq.(]7|) we drop the terms of higher order than 0(p) out. 
Thus, we come to the result for the 'kinematical' term k as 

32tt 2 

«=— . 12 

T 

But within the accuracy taken at calculating Eq. @ we are permitted to fix k at some point as 
k(p). 

Analysing the variations of the functional integration result when having calculated for 
the quasizero mode in the adiabatic approximation we see that because of the kinetic energy 
smallness it is certainly permitted to neglect its impact on the pre-exponential factor, which 
is small itself. The inverse infuence of the pre-exponential factor on the 'kinematical' term is 
negligible as well. Thus, going the way which we have already passed through while calculating 
Eq.(||) we receive whithin the required order of accuracy f] 

(S) = [d A z [ dp n(p) (~k p 2 + s(p)\ . (13) 



5 In principle, such an integral can be calculated exactly B. For the case at hand the types of deformations 
(quasizero modes) may be simply counted in terms of the instanton solution. They are induced by the variations 
of the instanton size, the changes of its position and colour orientation. In fact, there is one more type of the 
quasizero modes related to two far distant instantons but it has already been studied in H . 
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It is important to remark here that the contribution of PP pair interacting term averaged over 
the colour orientation will be 

(5(12)) = / ^ / ^ U M {12) , (14) 

tt n^_ 8f2 JL - to) 2 ] pM - ( \ -_i o 

^mH iZ J — — o~T72 7 / a X 4/ 2 , 2Y2 4? 2 , 2V2 ' Vi ~ X ~ Z i Pi ~ Pi\ z i) j I — 1, ^- 

# 2 iV c 2 - 1 7 2/i (2/i + Pi) 2 2/2(2/2 + P2) 2 

The way how one instanton of the size p% affects another of the size P2 is estimated by the 
following integral 

^U int (12) ~ p\{z 2 ) p 2 2 (z 2 ) . 

v g v 

Indeed, the adiabatic constraint makes it possible to rescale the integration variable ^ = 
d + Zpdpi ~ d (®* \ . Besides, we appraise the magnitude of a slowly varying function p\(z\) 
at the point where the integrand reaches its maximum. An analysis shows that the points z% 
and z 2 are not much separated as \z\ — z 2 \ ~ max{p 1 ,p 2 }, then the pi(z 2 ) looks like a quite 
suitable choice (the contact interaction) under the adiabaticity condition. In particular, if the 
PP sizes do not vary we reproduce the well-known result (Si nt (l2)) = y p\p\ for such a 

contribution. The integral / dp p 2 n(p)/n is approximated by p 2 in the approach we are going 
herein along (what is intuitively clear and can be easily argued) and eventually it results in the 
self-interaction of PP's. 

In the Minkowski space the factor in the curie brackets of Eq.(|l3|) might be interpreted as 
a mechanical system with Lagrangian C = | k p 2 — U e ff(p) and an action per one instanton 
might be taken as a 'potential energy' U e ff(p) = (3{p) + 51n(Ap) — ln/3 27Vc + v = . In the local 

vicinity of the potential minimum p 2 = ^ p 2 , (j-^f- = o) the system is oscillating and we 
have for the frequency (using the configuration corresponding to Eq. (|T2|) ) 

2 4z/ v 

K p p p z 



d?Ueff{p) I 4v 
dp 2 \pc p 2 ■ 



while calculating the second derivative 

We have only analysed the deformations in the temporary direction. Those in the spatial 
directions could be estimated by drawing the same arguments. Thus, the expression for the re- 
term keeps the form obtained above with the only change of rates for the appropriate gradients 
of function p(x), i.e. the substitution p{t) — > should be performed for such a 'crumpled' 
instanton. Then the frequency of proper fluctuations might be interpreted as the mass term 
and the excitations occur to have a phonon-like nature 

C= l -K[p 2 -V P Vp\-U eff {p) , (16) 

(the cross-terms ~ p p equal to zero identically) 0. The parameters p and [3{p) determined 
(by maximizing the partition function of II liquid with respect to N 0) in a self-consistent 

6 It is interesting to remark that the centre of instanton solution may not be shifted since the relevant 
deformations lead to the singular k (unlike the dilatational mode), whereas the colour coordinate variation Q 
gives the trivial result k = 0. 
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way take the following values pA rs 0.37, (3 ~ 17.5, n A 4 « 0.44, (iV c = 3) therefore, 
for the mass term we have m ~ 1.21A. Then the wave length A4 ~ m -1 in the a^-direction is 
A 4 A w 0.83 > y x / 4 A « 0.81 > pA (the size of the phonon localization in the spatial directions 
can be arbitrary and may noticeably exceed A4, besides the number N of the PP's forming the 
excitations might be pretty large). 

These numerical values obtained should be taken rather qualitatively illustrating the prin- 
ciple possibility to have the particle-like excitations originated by the quasizero modes. It is 
clear that striving to go beyond the superposition ansatz one should to take into account at 
least a medium change of the instanton profile and to develop more realistic description of the 
instanton interactions. Certainly, what presented here essentially exceeds the corresponding 
results which one could expect in the 'complete theory'. 

In conclusion let us emphasize that we have considered the excitations of the instanton 
liquid generated by the dilatational instanton deformations and the adiabaticity assumption 
leads, in principle, to a fully consistent picture. The model itself regulates the most suitable 
regime of such phonon-like deformations resulting in mass gap generation fixed by Aqcd- Ap- 
parently, including the quark condensate [] an intriguing guess is to associate some light hadrons 
with these phonon-like excitations discovered since the preliminary evaluations of their mass 
spectrum look quite encouraging. Moreover, the concept of the confining potential for the 
light quarks in the context of our approach seems simply irrelevant because of a stable phonon 
nature. 
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